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Abstract 

Given an anisotropic Voronoi diagram, we address the fundamental question of when its 
dual is embedded. We show that, by requiring only that the primal be orphan- free (have 
connected Voronoi regions), its dual is always guaranteed to be an embedded triangulation. 
Further, the primal diagram and its dual have properties that parallel those of ordinary Voronoi 
diagrams: the primal's vertices, edges, and faces are connected, and the dual triangulation has 
a simple, closed boundary. Additionally, if the underlying metric has bounded anisotropy (ratio 
of eigenvalues), the dual is guaranteed to triangulate the convex hull of the sites. These results 
apply to the duals of anisotropic Voronoi diagrams of any set of sites, so long as their Voronoi 
diagram is orphan- free. By combining this general result with existing conditions for obtaining 
orphan-free anisotropic Voronoi diagrams, a simple and natural condition for a set of sites to 
form an embedded anisotropic Delaunay triangulation follows. 



1 Introduction 



Voronoi diagrams and their dual Delaunay triangulations are fundamental constructions with nu- 
merous associated guarantees, and extensive application in practice [5l[7]. At their heart is the use 
of a distance between points, which in the original version is taken to be Euclidean. This suggests 
that, by considering generalizations of the Euclidean distance, it may be possible to obtain variants 
which can be well-suited to a wider range of applications. 

Attempts in this direction have been met with some success. Power diagrams [8j generalize 
Euclidean distance by associating a bias-term to each site. The duals of these diagrams are guar- 
anteed to be embedded triangulations, in any number of dimensions. Although this is a strict 
generalization of Euclidean distance, it is a somewhat limited one. The effect of the bias term is 
to locally enlarge or shrink the region associated to each site, loosely-speaking "equally in every 
direction". It allows some freedom in choosing local scale, with no preference for specific directions. 

Another way to generalize Voronoi diagrams is to endow Euclidean space with a continuously 
varying Riemannian metric, and use its associated (geodesic) distance in the definition of the 
Voronoi diagram. This generalization carries a significant amount of freedom, allowing for local 
scale and "directionality" to be freely specified at each point. 

Despite being of potentially great interest, this latter approach has faced several obstacles. The 
dual of such an anisotropic Voronoi diagram is, in general, not an embedded triangulation, even if 
the primal is orphan- free (the Voronoi regions are connected), and may produce element inversions 
and edge crossings. An additional obstacle is that the geodesic distance is extremely expensive to 
compute in practice, rendering traditional methods for constructing Voronoi diagrams impractical. 

Two notable attempts to resolve these problems have been proposed. Independently, Labelle 
and Shewchuk [12], and Du and Wang [6J propose two efficient approximations of geodesic distance 
from a site to any point in the domain. Although their associated Voronoi diagrams are, in general, 
no longer orphan-free, Labelle and Shewchuk show that a set of sites exists with an orphan-free 
diagram, whose dual is embedded, in two dimensions. They accomplish this by proposing an 
iterative site-insertion algorithm that, for any given metric, constructs one such set of sites. Some 
recent work making use of the above definitions are [21 [1] . 

In this paper, we show that if a set of sites produces an orphan-free anisotropic Voronoi diagram, 
using the definition of [6J , then its dual is always an embedded triangulation (or a polygonal mesh 
with convex faces in general), in two dimensions (Thm. 6.2). This effectively states that, regardless 
of the sites' positions, if the primal is well-behaved, then the dual is as well. Further, in a way 
that parallels the ordinary Delaunay case, the dual has no degenerate elements (Lem. 4.9), its 
elements (vertices, edges, faces) are unique (Cor. 6.3), and, under mild assumptions on the metric, 
is guaranteed to triangulate the convex hull of the sites (Thm. 4.7). Note that, while [12] prove a 
property of the output of an algorithm, the results in this paper are fundamental properties of the 
diagrams themselves (independent of how the orphan- freedom of the primal was obtained). 

In practice, we may combine our results with those of [3J to conclude that duals of anisotropic 
Voronoi diagrams of appropriate e-nets are always embedded triangulations (Cor. 6.4). This may be 
particularly useful in (asymptotically-optimal) function approximation applications, where we are 
often interested in constructing the anisotropic Delaunay triangulations of appropriate e-nets [4, llj, 
and was the initial motivation for the current work. Finally, we note that, as discussed in Sec. [7| 
in practice, algorithms for constructing anisotropic diagrams of the type discussed here may have 
a more numerical flavor (i.e. front-propagation, fast marching methods), in contrast with the more 
combinatorial nature of other Voronoi diagrams. 
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2 Setup 



Given a finite set V C M? oi sites on the plane, a Voronoi diagram decomposes Euclidean space 
into regions, each distinguished by the site its points are closest to. The notion of closeness is 
defined as follows: consider a continuous metric (in coordinates: Q : M? R^^^, symmetric, 
positive definite) over two-dimensional Euclidean space and, following [6J, define the (asymmetric) 
"distance" between a site {v G V) and point (p G M^) as 

Div,p)=[ip-vYQpip-v)Y^^ 
The Voronoi region of site v is the set of points no further from v than from any other site: 

R(v) = G : D{v,p) < D{w,p),yw G V} 

Note that, because y is a set, and not a multiset, there are no coincident sites. 

This definition of anisotropic Voronoi diagram results in diagrams that take the usual shape, 
but whose regions may have curved boundaries, and tend to be elongated along certain directions, 
depending on the metric Q. 

Following the notation of [12j, we say that a diagram is orphan-free if its regions R(v) are 
connected, for all v ^ V (equivalently: each cell, or connected component of a Voronoi region, 
contains its generating site). 

An important distinction in this construction is that the "interfaces" in-between Voronoi regions 
do not, in general, have null measure (and points equidistant to three or more sites are not always 
isolated), unlike geodesic-distance diagrams [13j (see ^ Lemma 5.2), and those of \12\ . 

Note that, given one such degenerate diagram, any interior point p of a Voronoi edge R{v)nR{w) 
of non-zero measure has an open neighborhood Np C R{v) H R{w). By adding to Q an arbitrarily 
small perturbation supported on Np, an orphan is created. That is, with respect to the metric, 
any degenerate diagram is arbitrarily close to a diagram that has orphans, and therefore the non- 
degeneracy requirement is only slightly more restrictive than the orphan-freedom one. 

We begin with three simple lemmas, proved in Appendix A, which provide an introduction to 
some of the techniques used in the sequel. The first is a simple extension of Lemma 2.1 of [6], and 
the second follows directly from the continuity of Q and the fact that every site is strictly closer to 
itself than to all other sites. 

Lemma 2.1. Given two sites v,w ^ V , v ^ w, the only point equidistant to v, w in their supporting 
line is the midpoint (v + w)/2. 

Lemma 2.2. Every site is an interior point of its corresponding Voronoi region. 

Lemma 2.3. Every Voronoi region of an orphan-free anisotropic Voronoi diagram in M? is simply 
connected. 

Primal diagram. The Voronoi diagram is a collection of Voronoi regions R{v), one for every site 
in V, as well as a structure induced by the sets of points closest and equidistant to two or more 
sites. Let P to be the primal Voronoi diagram embedded on the sphere by stereographically 
projecting the plane onto the punctured sphere and completing it with a "point at infinity" poo- P 
includes not just regions, but also embedded edg es Eynj, which are connected sets of points closest 
and equidistant to two sites v,w, as well as vertices (points closest and equidistant to three or 
more sites). Since we are assuming that sets of equidistant points to two sites have null measure, 
and that points equidistant to three or more sites are isolated, we can consider the abstract mesh 
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P — (V^, Ep^ Fp) with structure derived from P, where Vp, Ep, Fp are the vertices, edges, and faces 
of P. 

A simple induction argument reveals that every edge in Ep connects two vertices in V^. To see 
this, start with an orphan-free diagram with two sites. The only edge of E is the unbounded set 



of points equidistant to the sites, passing through their midpoint (Lem. 2.1), and with endpoints 
at Poo- The insertion of an additional site can only split existing edges of E at points where two 
edges cross (a vertex by definition), and thus the original edge is split into pieces whose endpoints 
are now either p^o or an edge crossing (a vertex). 

Clearly, P is a planar embedded diagram since two edges cannot meet except at vertices: any 
point p where edge D(u,p) = D{v,p) crosses edge D{v,p) = D(w,p) is equidistant to three sites, 
and thus a vertex. Therefore P is planar. 

Dual of an orphan-free Voronoi diagram. The dual of the primal mesh structure P is another 
abstract mesh G = (V,E,F), with one vertex per site, and edges connecting adjacent regions of 
P. Its face structure is derived from that of P by duality and, in particular, (V^ E) is planar, since 
(V^, Ep) is planar. Using the embedded primal P, we can define an embedding G of G on the plane, 
with vertices coinciding with sites (and curved edges). Although P was defined on the sphere, we 
embed G on the plane instead, and therefore we exclude (from both G and G) the unbounded face 
dual to the point at infinity Poo- 

Since (y, E) is planar, G derives the following properties from the primal: 

(i) Every edge in E is incident to two dual faces, except for boundary edges (whose corresponding 
primal edge in P connects unbounded regions), which are only incident to one dual face. 

(ii) To every dual face in P, connecting vertices vi^ . . . ^Vm^ corresponds a primal vertex c G 
R{vi) n • • • n R{vm) of P, which is equidistant to the corresponding sites . . . , Vm. 

The property (ii) requires that the Voronoi diagram be orphan-free, as noted in the remark after 
Thm.O 

Notice that either G or P could be multigraphs (e.g. P would be a multigraph if an edge Ep 
of P has multiple connected components). However, the following lemma ensures that both are 
simple graphs. 

Lemma 2.4. G and P are simple (have no multi-edges or self-loops). 

Proof. See Appendix D. □ 

By its definition, an edge of P is a connected set of points equidistant and closest to two sites 
v^w. Since P has no multi-edges, this implies the following: 

Corollary 2.5. The set of points closest and equidistant to two sites of an orphan-free diagram is 
connected. 

Property (ii) above can be restated in a more useful way that parallels an equivalent property 
of duals of ordinary Voronoi diagrams. 

Definition 2.6 (Empty circum-ellipse property). A face of G incident to vertices vi,. . . ,Vm 
satisfies the empty circum-ellipse (ECE) property if there is some ellipse that circumscribes the 
sites corresponding to vi, . . . ,Vm, and contains no site in its interior. 

Note that any empty circum-ellipse serves as witness to the ECE property. We show that the 
ellipse centered at a Voronoi vertex c, with axis given by Qc, is always an empty circum-ellipse of 
its dual polygon. 
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Theorem 2.7. Every face of G satisfies the empty circum- ellipse property. 

Proof The proof is constructive. To every face in G connecting vertices with corresponding sites 
vi,. . . ,Vm, by property (ii), corresponds a Voronoi vertex c G R{vi) H • • • H R{vm) of P that is 
equidistant to vi^ . . . ^Vj^- If /i = c) = • • • = D{vm^ c), then the open ehipse 

^ = e : D{p, c) < = E : - c)^Qc{p - c) < m^} 

circumscribes the vi^. . . and does not contain any site (otherwise, c would be closer to that 
site than to . . . , and therefore not in R{vi) H • • • H R{vm)^ a contradiction). □ 



(a) 



(b) 



Figure 1: In an orphan- free diagram (a), every face satisfies the empty-circum ellipse (ECE) prop- 
erty. If the diagram is not orphan-free, then some faces may not satisfy the ECE property. 

Although this may not at first seem apparent, the above property uses the fact that the primal 
diagram is orphan free. The critical fact in the proof of Thm. |2.7| is not that there is a Voronoi 
vertex c equidistant to i^i, . . . , but rather that it is c G R{vi) H • • • H R{vm). Consider Figs. 1(a) 
and |l(b){ in which sites vi^ vj^ Vk form a face of G. In 1(a) , the diagram is orphan- free, and the fact 
that Cijk ^ R{vi) n R{vj) n R{vk) implies Thm. |2.7[ Figure 1(b) shows part of a Voronoi diagram 
that has an orphan cell (corresponding to site w) "covering" Cijk- If the dual is constructed such 
that the face Vi,Vj,Vk is not in the dual then a "hole" in the triangulation may occur. If the face 
Vi^ Vj^Vk is in the dual, then, despite the fact that Cijk is equidistant to Vi^ vj^ v^^ this Cijk is no longer 
in R{vi) n R{vj) n R{vk)^ but in the Voronoi region of w. Since cijk is closer to w than to Vi^ vj, Vk, 
any open ellipse centered at Cijk, and having Vi^vj^Vk on its boundary, contains w. Therefore, in 
Fig. 1(b), the ECE property does not hold for the dual face of Cijk- Ensuring that Thm. 2.7 holds 
is one of the main reasons for requiring orphan-freedom. 



3 Summary of Results and Outline 

The two main results in this paper assume that we are given an orphan-free anisotropic Voronoi 



diagram (in two dimensions) and prove that: 1) the dual is an embedded triangulation (Thm. 6.2) 



and 2) if the metric has bounded anisotropy (ratio of eigenvalues), the dual triangulates the convex 



hull of the sites (Thm. 6.1). 



We prove a number of additional results. Possibly the most important shows that the elements 
of an orphan- free diagram (vertices, edges, faces), are connected sets, and therefore, in some sense. 



unique (Cor. 6.3). This is important, since this is a natural well-behave-ness condition on the primal, 
which, together with the above results, indicates that orphan- freedom is sufficient to guarantee well- 
behave-ness of both primal and dual. 
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A smaller result of interest is: the dual faces satisfy an empty cir cum- ellipse property (Tlim.[2^, 
which parallels the empty-circumcircle property of ordinary Voronoi diagrams, and could have 
further practical implications from the ones described here. 

The order of the proofs is, however, different from the one just stated. In effect, we begin by 
proving the most restrictive case: that the dual of an orphan-free diagram with metric of bounded 
anisotropy is embedded (Sec. [5]), and that its boundary is the boundary of the convex hull of the 
sites (Sec. [i]). We then relax the bounded anisotropy condition (Sec. |6]), and show that, in the 
general case, we loose the convex hull property, but the dual remains embedded. Finally, once 
these results are established, we prove that vertices of the primal are unique (uniqueness of edges 
was proven in Sec. [2) and uniqueness of primal faces is the same as orphan- freedom) . 

We assume in the remainder of the paper that there are more than two sites, not all of which 
are colinear, and relegate the (considerably simpler) colinear case to Appendix B. 

4 Dual of Orphan-free Diagram (Part I: boundary) 

In this section, we assume that the metric has bounded anisotropy, and conclude that the boundary 
of the dual of an orphan-free diagram is the same as the boundary of the convex hull of the sites 
(and in particular is simple and closed). If RpApR^ is the eigendecomposition of (5^9 at p G M^, 
with Ap = diag[X2{p)^ Xi{p)], X2{p) > M{p) > 0, and Rp orthonormal, then we assume that there 
is some bound 7 on the anisotropy of Q, such that: 1 < X2{p) / Xi{p) < 7^ for p G M^. Note that 
this condition may commonly hold in practice, if the metric is sampled on a compact domain (and 
possibly extended to the plane by reusing sampled values only). 

We begin by defining G to be the straight-edge drawing of G with vertices at the sites. For 
the moment, we assume that every Voronoi vertex is equidistant to no more than three sites, and 
therefore that all faces in G are triangles (we extend the results to the general case in Appendix G). 
We associate to (5 a mapping (from the vertices, edges, and faces of G to M?) for which, because 
all faces of G are triangles, it is well defined whether a point in M? belongs to any given face of 
(5, a fact that will be used in the proofs of Sec. [sj This mapping will be shown in Sec. |5] to be 
an embedding. In the sequel, it is assumed that G encompasses both the mesh structure, and the 
mapping. 

The boundary vertices of G are those whose corresponding primal regions in P are unbounded, 
while boundary edges of G connect boundary vertices. Note that this is a topological property of 
G, rather than a geometric one (boundary elements of G may, in principle, not lie in the boundary 
of the convex hull of the sites). For convenience, we call B C E the set of boundary edges of G. 

The convex hull CT-L(V) of V is the minimal (w.r.t. set containment) convex set that contains 
V. For convenience, we name W = {wi G : i = 1, . . . , m} the sites that are part of the boundary 
of the convex hull CT-L{V), and order them in clock- wise order around CT-L{V). The boundary B 
of the convex hull is a simple circular chain B — {{wi^Wi^i) : z = 1, . . . ,m}. We prove that it 
IS B — B (loosely speaking: the topological boundary of G, and the geometric boundary of its 
straight-edge embedding G, are the same), which implies that G covers the convex hull of the sites, 
and its boundary edges form a simple, closed polygonal chain. All the proofs of this section are in 
Appendix E. 

Lemma 4.1. To every boundary edge (v,w) of G corresponds a segment in the boundary ofCT-LiV). 
[B(ZB] 

We now turn to the converse claim: that to every segment {wi^ wj) ^ B corresponds a boundary 
edge (wi^Wj) G S in G. Since B is the set of boundary edges of G, whose primal edges in P are 
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unbounded, the claim is equivalent to proving that, to every segment (wi^wj) in the boundary of 
the convex hull corresponds an edge of the dual G, whose primal edge E^.^^^j is unbounded. 

The proof proceeds as follows. First, assume w.l.o.g. that the origin is in the interior of CT-L{V). 
Let C{a) = {x G : ||x|| = a} be a sufficiently large origin- centered circle. We define two 
functions: 

7T : C{a) dCH{V)^ 7t{p) = argmin D{r^p) 

recn{v) 

V : dCU{y) z/(r) = a • r/||r|| 

TT simply projects every point in C to its closest in Cl-iiy^ (with respect to the distance I?), and v 
projects a point back to C. 

Lemma 4.2. — argmin is a continuous function in 'E? . 
xecn{v) 

Proof. We first prove by contradiction that 7v{p) is unique. Let Mp be the unique symmetric 
positive definite square-root matrix such that M^M^ = Qp. Consider distinct ^1,^2 ^ CHiV) 
closest to p. Since ^1 ^ ^2, by the convexity of the Euclidean norm, and the positive-definiteness 
of M^, it is 

D{{qi + q2)/2,p) = \\Mp {{qi + q2)/2 - p)) \\ < \\Mp {qi - p) \\ = \\Mp - p) \\ 

Thus p is closer to {qi + ^2)/2 than to qi, q2, a contradiction. Therefore the closest point 7v(p) to p 
is unique. 

If TV were not continuous at p, then there is e > such that for all 5 > there is a point p = p{6) 
such that \\p — p\\ < S and \\7t{p) — 7t{p)\\ > e. Consider the sequence {^i = 7r(p(l/i)) : z G N} 
of points in CT-L{V). Because CT-L(V) is compact, {nci} has a subsequence that converges to some 
G CTi{V). By continuity of D (which follows from the continuity of Q), it is 

D{7v{p),p) = min D{x,p) = lim min D(x,pi) = lim D{Hii,pi) = D{hc,p) 

xeCH(V) i^ooxeCH(V) i^oo 

and therefore the closest point in CT-L(V) to p is not unique, a contradiction. □ 

Clearly, u is continuous in dCT-L{V). Note that, because CT-L{V) contains the origin, then, as 
shown in Fig. [2j u projects every point 7t{p) G {wi^Wj) on a segment of dCT-L{V), outwards from 
the convex hull (and on the empty side of (wi^Wj)); that is, so that u{7t{p)) G D C (z/(7r(p)) is 
in the empty half-space of ( wj, Wj)). 



Conveniently, Lem. 4.3 shows that, if 7r{p) = {wi + Wj)/2^ then p G E^.^^.^ and so the claim 
now reduces to showing that for each segment {wi,Wj) of dCT-L{V), and for every sufficiently large 
circle C, there is j9 G C with T^{p) — (tt^i H" Wj^jl. Since this implies that E^^^^. is unbounded, it 
means that the corresponding edge {wi,Wj) is in 

The proof is by contradiction. Lem. |4.5| uses Brouwer's fixed point theorem to show that, for 
every segment (wi^wj) of B, if there were no p G C with 7v{p) = {wi + Wj)/2^ then the function 
V o 71 \ C ^ C must have a point q ^ C such that u{7r(q)) = — or equivalently such that q is 
"behind" the segment (wi^wj) G dCT-L{V) to which it is closest {q G H^). On the other hand, 
Lem. |4.4| shows that, for all sufficiently large circles C, no point g G C can be closest to some 
segment {wi,Wj) G dCT-L{V) that it is behind of, creating a contradiction. 

As mentioned above, we are interested in identifying points in C{a) that are also in some primal 
edge Er^.^^.. The following lemma will be used to show that, for sufficiently large a, such a point 
p can be characterized by satisfying 7r(p) = (wi + Wj)/2. 
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Figure 2: The construction for the proof of Thm. 



4.7 



Lemma 4.3. If 7v{p) = (wi ^Wj)/2, with p G , Wi,Wj G W, then D{wi,p) = D(wj,p). 

The following technical lemma is the key in constructing a contradiction by showing that, for 
sufficiently large circles C(cr), no point q G C{a) can be closest to some segment {wi, wj) E dCT-L(V) 
that it is behind of {q ^ Hij)i where, as before, the open half space is chosen to be the only of 
the two half spaces on either side of the supporting line lij of (wi^ wj) such that nV ^ (j). 

Lemma 4.4. There is p such that, for any segment (wi^wj) G B, with supporting line lij, every 
p G with \\p\\ > p whose closest point in kj is rrip G WjUJ] is closer to a site inV\ {wi, wj} than 

to lij. 

Lemma 4.5. Every continuous function F : ^ that is not onto has a fixed point. 

Lemma 4.6. To every segment in the boundary of CT-LiV) corresponds a boundary edge of G. 
[B D B] 

Therefore, 

Theorem 4.7. If the metric Q has bounded ratio of eigenvalues, then the boundary edges of G are 
the same as the boundary edges of the convex hull ofV. [B — B] 

Corollary 4.8. // the metric Q has bounded ratio of eigenvalues, then the boundary edges of G 
form a simple, closed, convex polygonal chain. 

The ECE property is sufficient to show that no face in G is degenerate, a fact that will be used 
in the following section. 

Lemma 4.9. G has no degenerate (null area) elements. 

In Thm. |5.9[ we show that, even if Voronoi vertices are incident to more than three Voronoi 
regions (equidistant to more than three sites), every face in (5 is a (strictly) convex polygon, and 
therefore can be trivially triangulated (say in a fan arrangement). The resulting triangulation has 
no degenerate elements, and all its triangles satisfy the ECE condition (using the same witness 
ellipsoid as the convex polygon from which they are triangulated). 
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5 Dual of Orphan-free Diagram (Part II: interior) 



In this section, we assume that the (topological) boundary of G is simple and closed, and prove 



that G must be embedded. The main argument in the proof uses Theorems |2.7| and |4.7[ Cor. |4.8 
as well as the theory of discrete one- forms on graphs, to show that there are no "edge foldovers" 
in G (edges whose two incident faces are on the same side of its supporting line), and use this to 



conclude that G is embedded (Thm. 5.9). As in Sec. [4| we assume that not all sites are colinear 



(the simpler cohnear case was addressed in Appendix B). We distinguish between the sites W 
that lie on the boundary of the convex hull, and the remaining, or interior sites {V \ W). 

The following definition, from [lOj, assumes that, for each edge (i;^, Vj) in G, we distinguish the 
two opposing half-edges {vi^Vj) and (vj^Vi). 

Definition 5.1 (Gortler et al. [lOj). A non-vanishing (discrete) one-form ^ is an assignment of a 
real value ^ij ^ to each half edge {vi^ vj) in G, such that ^ji = —^ij- 

Since G has the same structure as G, we can construct a non-vanishing one-form over G as 
follows. Given some unit direction vector n ^ S>^ (in coordinates n = [^1,712] ), we assign a real 
value z{v) = n^v to each vertex v in (5, and define ^ij = z(vi) — z(vj), which clearly satisfies 
^ji = —^ij- The one-form, denoted by is non- vanishing if, for all edges (vi^vj) G it is 
^ij = n}{vi — Vj) 7^ 0. That is, if the direction n is not orthogonal to any edge. The set of edges 
has cardinality \E\ < \V\{\V\ — l)/2, and in particular it is finite. Therefore almost all directions 
n G S-^ generate a non- vanishing one- form 

Since G = {V, F) is an planar graph with a well-defined face structure, there is, for each face 
f ^ F, Si cyclically ordered set df of half-edges round the face. Likewise, for each vertex v G V, 
the set 6v of cyclically ordered (oriented) half-edges emanating from each vertex is well-defined. 

Definition 5.2 (Gortler et al. ^IQJ). Given non-vanishing one-form corresponding to n ^ , 
the index of vertex v with respect to is ind^^(i;) = 1 — sc^^('u)/2^ where sc^ri{v) is the number 
of sign changes of ^'^ as one visits the half-edges of Sv in order. 

The index of face f is ind^^(/) = 1 — sc^^(/)/2 where sc^^(/) is the number of sign changes of 
as one visits the half-edges of df in order. 

Note that, by definition, it is always ind^^(i;) < 1. A discrete analog of the Poincare-Hopf 
index theorem relates the two indices above: 

Theorem 5.3 (Gortler et al. [lOj). For any non-vanishing one-form , it is 

^ind^.(i;)+^ind^4/) = 2 

v^v feF 

Note that this follows from Theorem 3.5 of [10] because the unbounded, outside face, which is 
not in G, is assumed in this section to be closed and simple, and therefore would have null index. 
Note that the machinery from [10] to deal with degenerate cases isn't needed here because vertices, 
by definition, cannot coincide (V is not a multiset). All proofs in this section, except for that of 



Thm. 5.9, are in Appendix F. 



The one-forms defined above satisfy the following property. 

Lemma 5.4. Given a non-vanishing the sum of indices of interior vertices (V \ W) of G is 
non-negative. 

The next two lemmas relate the presence of edge foldovers and the ECE property of Defini- 



tion 2.6 to the indices of vertices in G. 
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Lemma 5.5. // G has an edge foldover, then there is a non-vanishing one-form such that 
ind^n (v) < for some interior vertex v \ W . 

Lemma 5.6. Given n G S-^ and non-vanishing one-form , if G has an interior vertex v ^V\W 
with index ind^^(i;) = 1, then there is a face f of G that does not satisfy the empty cir cum- ellipse 
property. 

The above provides the necessary tools to prove the fohowing key lemma. 

Lemma 5.7. G has no edge foldovers. 

Finally, the absence of edge foldovers, together with a simple and closed boundary, is sufficient 
to show that G is embedded. 

Lemma 5.8. If its (topological) boundary is simple and closed, then the straight-line dual G of an 
orphan-free diagram, with vertices incident to at most three sites, is an embedded triangulation. 

As shown in Appendix G, even if there are non- generic Voronoi vertices that are incident to 
more than three sites, the dual is composed of faces each of which is convex and satisfies the ECE 
condition (Def. |2.6[ ). Each convex face can be triangulated (e.g. in a fan arrangement) in such a 
way that individual triangles satisfy the ECE condition with the same witness ellipse as the face 
from which they are derived. This leads to the following: 

Theorem 5.9. If its (topological) boundary is simple and closed, then the straight-line dual G of 
an orphan-free diagram is an embedded polygonal mesh with convex faces. 



6 Final Results 



We can combine Corollary |4.8| and Theorem |5.9| into 

Theorem 6.1. // the metric Q has bounded ratio of eigenvalues, then the dual of an orphan-free 
Voronoi diagram with respect to Q is an embedded polygonal mesh with convex faces, and covers 
the convex hull of the sites. 

This result can be generalized, dropping the bounded anisotropy condition on Q, but at the 
cost of losing the convex hull property, as shown in Appendix H: 

Theorem 6.2. The dual of an orphan-free Voronoi diagram is an embedded polygonal mesh with 
convex faces. 



From the above results. Corollary 2.5, and the definition of the dual, it follows (see Appendix 
H) that 

Corollary 6.3. An orphan-free anisotropic Voronoi diagram is composed of unique (connected) 
vertices, edges, and faces. 



Finally, we note that Thm. 5.9 can be combined with existing conditions for orphan-freedom [3], 
resulting in a simple and natural condition for a set of sites to induce an embedded polygonal mesh 
as the dual of their anisotropic Voronoi diagram: 

Corollary 6.4. IfV is an asymmetric e-net w.r.t. D, Q a continuous metric with metric variation 
a, andea < 0.09868^ then the dual of the anisotropic Voronoi diagram ofV is an embedded polygonal 
mesh with convex faces. 



10 



where an asymmetric e-net is simply a weaker form of e-net defined on non-symmetric functions 
which can be computed with the iterative algorithm of [9J, and the metric variation a is a 
Lipschitz-type condition on Q [3]. The above condition is known to be conservative, and there may 
be simpler conditions to achieve orphan-freedom. As a practical observation, Du and Wang [H] 
report orphans to be a rare occurrence in their experiments. 



7 Proof-of-concept Implementation 




(a) 



(b) 



(c) 



(d) 




(e) 



(f) 



(g) 



(h) 



Figure 3: Anisotropic Voronoi diagrams, and their duals generated by our proof-of-concept imple- 
mentation. Voronoi vertices are marked in red, while dual vertices (sites) and edges are drawn in 
black. 

Though not aiming for an efficient implementation, we implemented a simple proof-of-concept 
that constructs anisotropic Voronoi diagrams like the ones considered in this paper, and their 
duals (Fig. 3). A closed- form metric, which has bounded ratio of eigenvalues, is discretized on a 
fine regular grid, and linearly interpolated inside grid elements, resulting in a continuous metric. 
The sites are generated randomly (Figs. |3(a)| and |3(b)[ ), or using a combination of random, and 
equispaced points forming an asymmetric e-net (remaining figures). 

The primal diagram was obtained using front propagation from the sites outwards, until fronts 
meet at Voronoi edges. The runtime is proportional to the grid size, since every grid- vertex is 
visited exactly six times (equal to their valence). The implementation does not guarantee the 
correctness of the diagram unless it is orphan-free, and serves to verify the claims of the paper 
since well-behave-ness of the dual is predicated on that of the primal. 

The two main claims of the paper are clearly illustrated in these examples. In all examples, the 



dual covers the convex hull of the vertices (Thm. 4.7), is a single cover, embedded with straight 



edges without edge crossings (Thm. |5.9| ), and has no degenerate faces (Lem. 4.9). By focusing 
on the primal diagrams (second and fourth column), further claims in the paper become apparent. 



namely that Voronoi regions are simply connected (Lem. 2.3), Voronoi vertices are unique (Cor. 6.3) 



Voronoi edges are connected (Cor. 2.5), unbounded regions correspond to boundary dual vertices, 
and unbounded edges of the Voronoi diagram correspond to boundary dual edges. 
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8 Conclusion and Open Questions 



We studied the properties of duals of orphan-free anisotropic Voronoi diagrams, for the purposes of 
constructing triangulations on the plane. The main result (Theorems 5.9 and 4.7, Cor. 6.4) is that 
the dual, with straight edges and having the sites as vertices, is embedded and covers the convex 
hull of the sites, mirroring similar results for ordinary Voronoi diagrams and their duals. 

A few, somewhat less important properties are proven, including the fact that every primal 
region is simply connected, that elements of the primal are unique (Cor. 6.3, but perhaps most 
significantly that every face in the dual satisfies an empty circum- ellipse property that has a direct 
parallel in the empty circum-circle property of ordinary diagrams, and is the basis for proving that 
it is embedded with straight edges. 

Perhaps the most important outstanding question may be whether these ideas extend to higher 
dimensions. The results in Sees. |4] and [5j except for Lem. 2^, can be trivially extended to n 
dimensions. Sec. |4] has been written only for the two-dimensional case, but a similar construction, 
and the same argument would work in higher dimensions (Lem. 4.5 being a hint of this). It is the 
argument in Sec. [5] that becomes problematic. While the ECE property is shown to be sufficient 
to prevent foldovers in the triangulation, it is not sufficient in higher dimensions. In particular, 
fixing the boundary to be simple and convex, there are simple arrangements of tetrahedra in 
that contain face foldovers but do not break the ECE property. We plan to study these question 
next. 
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Appendix A 




Figure 4: Diagram for the proof of Lem. 2.3 



Lemma |2.1| Given two sites v,wgV,v^w, the only point equidistant to w in their supporting 
line is the midpoint (v -\- w)/2. 

Proof. If p is in the supporting hne of then p = (1 — X)v + Xw, A G M. p is equidistant to w 
(i.e. D{v,p) = D{w,p)) if and only if \^{v - wYQp^v - w) = {1 - X)'^{v - wyQp{v - w). 

Since Q is positive definite and v it is {v — wyQp(v — w) > 0^ and so the above equahty 

holds iff A = 1/2 is the midpoint). □ 



Lemma |2.3| Every Voronoi region of an orphan-free anisotropic Voronoi diagram in R is simply 
connected. 

Proof. A multiply connected region R{v) C M? is path connected, and is such that there is a map 
/ : ^ R{v) that cannot be continuously contracted to a point. We can assume / injective 
(simple), since a non-injective / can always be broken up into injective pieces, at least one of which 
must be such that it cannot be continuously contracted to a point (or else / would be). By the 
Jordan curve theorem, / encloses a bounded set B. Since / cannot be continuously contracted to 
a point, there must he I C B with / C \ R{v). Since / is a subset of B, and B is bounded, / is 
bounded. 

/ is part of the Voronoi diagram, so it must be composed of one or more Voronoi cells. Because 
the diagram is orphan-free, each cell in the Voronoi diagram contains its generating site: / = 
UiR{wi)^i = 1, . . . , m (Fig |4(a)'] ), for some number m > 0. 

We can now consider what the diagram would look like if we remove, one by one, all of the 
sites Wi in /, until only one (wi) is left. The new set of sites is ^ = y \ {wi : i = 2, . . . , m}. From 
the definition of Voronoi diagram it is clear that the region corresponding to site wi in the new 
diagram is R'{wi) C / (since no point is strictly closer to wi than to {wi : i = 1, . . . ,m}), and 
R'{wi) 7^ (j) (since, by Lem. |2.2[ it is wi G int R'{wi)). 

In particular, since R'{wi) ^ 0, and R'{wi) is connected, R'{v) is still multiply-connected, with 
R'{wi) playing the role of / in the new diagram (Fig. 4(b)| ). We show that this cannot be the case, 
resulting in a contradiction. 

Since / is bounded, then R'{wi) C / must be bounded as well. And since wi is an interior point 
of R'{wi)^ and v is an interior point of R'{y) (and therefore an interior point of the complement of 
R'{wi))^ then the line passing through v and wi must intersect R'{v) at least at two distinct points 
p and q. 

Since the boundary between R'{v) and R'{wi) is composed of equidistant points to v^wi^ then 
both p^q must be equidistant to v^wi. By Lem. |2.1| only one point on the line connecting v^wi 
can be equidistant to them, a contradiction. 

□ 
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Appendix B: Colinear sites 




Figure 5: If all sites are colinear, the dual is always a polygonal chain. 

If all the sites are colinear then the structure of the Voronoi diagram P is greatly simplified, and 
always has the form shown in Fig. 5. In particular, P has no vertices aside from p^o since vertices 
are equidistant to three or more sites, and no point p E can be equidistant to three colinear 
sites (since points equidistant to p form an ellipse, and a line intersects an ellipse at most twice). 

Consider the set of sites V — {wi = ,m} ordered linearly along their supporting line. 
We shown that the graph dual to P, having V as vertices, has edges {{wi^ wij^i) = ..,m — 1}. 

For every pair of sites wi,wi^i, because they are consecutive, and by Lem. |2.1[ the midpoint 
mi^l^i = (wi + wi^i)/2 is closest and equidistant to wi^wi^i^ and therefore the edge (wi^wi^i) is 
in the dual. 

We finally show that every dual edge (wi^wj) is of the form (wi^wi^i). Assume otherwise, and 
therefore that, since Wi^wj are not consecutive, there is some between them. Because there is 
a dual edge (wi^wj), the corresponding primal edge E^.^^. in P is not empty, and therefore there 
is some point pij G Ey^.^y^. that is closest and equidistant to Wi^wj. By the convexity of D{-^pij)^ 
and the fact that Wk G WjUJ], it is D(wk^Pij) < D(wi^pij) = D(wj^pij)^ a contradiction. Therefore 
{wi,Wj) are consecutive sites. 

Since we have characterized the set of edges and vertices, and there are no (interior) faces in 
the dual, this completely determines the dual when all sites are colinear. 

Perhaps interestingly, this means that, in the colinear case, the structure of the dual does not 
depend on the metric. 

Appendix C: Technical Lemmas 

We include in this appendix all claims that, although needed to prove the results in this paper, 
are used only as intermediate steps. Some of them reveal useful aspects of the structure of the 
problem. In particular, Lem. |6.4] is used as a basic step throughout the paper. The proofs are quite 
technical, and can be skipped on a first read without affecting the comprehension of the remainder 
of the paper. 

The following result states that, even though the regions and edges of an anisotropic Voronoi 
diagram may be unbounded, there is a sufficiently large ball outside which there are no vertices. 

Lemma 6.1. The set of vertices in the primal diagram P is bounded. 

1 /2 

Proof. Let Mp = RpAp Rp be the unique symmetric, positive definite square root of Qp. Define 
Qp — Qp ' ^i^{p) — ^pdiag [A2(p)/Ai(p), 1] Rp (which is continuous by virtue of the continuity of 
Q and Xi{p)) and its symmetric positive definite square root Mp. 
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Pick three sites Vi, vj, G V . If they are cohnear then no point is equidistant to them, and so 
the lemma holds vacuously. If they are not colinear, consider some Q' G M^^^ symmetric, positive 
definite with ratio of eigenvalues in the range [1,7^]. Define ^^(a, 6, c) to be the unique center of 
a circle passing through three (non-colinear) points a, 6, c. ^ is continuous wherever a, 6, c are not 
colinear. 

If p is equidistant to Vi^ Vj,Vk^ and has = Q\ then clearly 

{vi - pfQ'ivi -p) = {vj - pYQ^Vj -p) = {vk - p)^Q\vk - p) 

or equivalently \\M'vi - M'p\\ \\M'vj - M'p\\ \\M'vk M'p\\. This means that M'p 
-^{M'vi, M'vj.M'vk), and so p = M'~^'^{M'vi, M'vj.M'vk) is the unique point that can be equidis- 
tant to Vi^ Vj, Vk while having Q'^ — Q'. 

Since the space of Q' is compact and M'~'^^[M' Vi^M' v^^M' v^) is continuous w.r.t. M' , then 
set of possible equidistant points to vi, vj, Vk is compact, and therefore bounded. Since this is true 
for every triple Vi, vj, Vk G V , and there is a finite number of triples, the set of possible equidistant 
points to three or more sites is compact, and therefore bounded. Since, by definition, vertices are 
equidistant to three or more sites, the lemma follows. □ 

For the sake of conciseness, we define: 

Definition 6.2. The open ellipse centered at p with m in its boundary is 

ep{m) = {x G : D{x,p) < D{m,p)} = {x G : {x-pYQpix-p) < {m - pYQpim - p)} 

Given some p ^ M?, we make frequent use of a map that transforms point s q G M? into 

We follow 



6.1 



= M^q, where = i?p diag[(A2(p)/Ai(p))-^/^, l]Rp defined in the proof of Lem 
the convention of denoting by q^ the transformed version of point q, where it will be clear from 
the context which matrix is used in the transformation. Since the eigenvalues of Mp are in the 
range [1,7], it's easy to see that: 

• Given p,q ^ M?, it is \\p — q\\ < \\p^ — Q^\\ ^ jWp — q\\ 

• // d(p,l) > is the minimum Euclidean distance between a point p and a line I, then 
7-'<d{p\n/d{pJ)<j. 

Proof: Given a,6 G /, a ^ 6, it is d{p',r) = \{p' - h') x {a' - h')\/\\p' - h'\\ ■ \\a' - h'\\ where 
the numerator is \{p' — h') x [a' — y)\ = 7|(p — x (a — and the denominator satisfies 
\\p-b\\-\\a-b\\ < \\p' -b'\\'\\a' -b'W <-i^\\p-b\\'\\a-b\\. Therefore /^Mp, ^ [7"^7]• 



The following technical lemma is the basis for the proof of Lem. 6.4 



Lemma 6.3. Given v,vu ^ M?, with supporting line I, and an open half space on one side of 
I, as well as a set S C of points whose closest point in I belongs to vw, then for all q G , if 
p ^ S has 

7^ max \\q — m\\ 
\\p\\ >max{||i;||, 11^11}+ '^^'''^ 



2d{qJ) 

then p is strictly closer to q than to I. 
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(a) (b) 
Figure 6: Diagram for the proof of Lem. |6.3 



Proof. Given a point p that satisfies the above constraint, by the definition of 5, it is rrip = 
argmin D{m,p) G vw. Also, since it is g G and since doesn't include q ^ I and therefore 

d{qj) > 0. Because rup = argmin D(m,p) G vw, it can be easily verified that the ellipse 9p{mp) 

must be tangent to / at rrip (Fig. 6(a)). When transforming it by M^, 9p{mp) becomes the circle 



'pi{m'p) = {x' G : \\x' —p'\\ < ||m^ ^11} (Fig. 6(b)), which is tangent to at m^. The lemma 



p' ^ . 

then reduces to showing that G 9'^,{m'^) since this implies q G Op(mp) and this in turn means that 
D(q,p) < D{mp,p) — min^^/ D{m,p), and therefore q is closer to p than to /. 

A simple calculation shows that q' G 0'^,{m'rp) iff ||y — m^|| > — m^||/ [2d{q' , I')]. Since rUp G vw 
implies \\mp\\ > max{||i;||, it is 

y - ^p\\ > \\P- ^p\\ > \\P\\ - W^pW > 
7^ max \\q — m\\ n n ii^/ _ ^/ ii 

mevw y ^2 11^ ^pll y 11^ ''^pll 



2d{q,l) - ' 2d{q,l) " 2d{q\U) 

and so q^ G 6p/{mp), as claimed. □ 

Given 1;^,^;^ G V whose Voronoi regions are neighbors in P, define and to be the two 
open half spaces on either side of the supporting line kj of Vi^Vj. Clearly, {H^ , ,lij} form a 
partition of M?. If the edge Ey.^y. of P corresponding to sites Vi,Vj is unbounded then, it must be 
either Ey.^y.nkj, E^.^y.^Hf-, or Ey.^y.nH^j unbounded. By Lem. 2.1, it is Ey.^y^nkj C {(vi + Vj)/2} 



and thus Ey.^y. nkj is always bounded. Therefore it must be that either Ey.^y. HH^ or Ey.^y. HH^ 
are unbounded (or both). The following lemma proves that to each unbounded Voronoi edge 
corresponds at least one open half-space that does not contain any site. 

Lemma 6.4. Given an edge Ey.^y. of the Voronoi diagram P, corresponding to neighboring sites 
Vi, Vj G V, if Ey.^y. n (Evi,vj n H^) is unbounded, then it is DV = (j) (H^ DV — (j)), where 
H^^H^ are open half spaces on either side of the supporting line ofvi^Vj. 



Proof Assume w G D V. We make use of Lem. |6.3| with Vi, vj , kj , Ey.^y. D taking the role of 
v,w,l,S, respectively. S ince E y^^y- OH^ is unbounded, we can choose p G Ey^^y. OH^ of sufficiently 
large norm, according to Lem. |6.3| such that D{w,p) < D{{wi + Wj)/2,p) which, by the convexity 
of D{',p), means that D{w,p) < D{{vi^Vj)/2,p) < D(vi,p) = D(vj,p). Since p is closer to w than 
to Vi, Vj, it is p ^ Ey.^y., a contradiction. □ 

Lastly, we show that every point of sufficiently large norm is strictly closer to the sites W = 
V n dCT-l{V) in the boundary of the convex hull of V, than to sites {V \ W) in its interior. This 
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(a) 



(b) 



(c) 



Figure 7: Diagram for the proof of Lem. |6.5 



means that, outside of a sufficiently large ball, the Voronoi diagram of V is the same as that of W, 
a fact critical in proving the results of Sec. [Ij The proof of this property makes use of the bound 
in the ratio of eigenvalues of Q, and it can easily be shown that counter-examples exist in cases 
where no such bound exists. The proof assumes that the origin of M? is chosen to be in the interior 

of cn{v). 



Lemma 6.5. There is p such that all p G 



with 



> p are closer to W than to V \ W . 



Proof. Begin by choosing p > max-^^y \\v\\ so that, since the origin is in int CT-L{v) 
norm greater than p is outside C1-L{V). Pick p such that it is also: 



every point of 



p > max \\w\\ + max 

wew (wr,ws)eB,wtew 



Y max^e^^j-^ \\x-wt\ 



2d{wt,l 



p > max I 

vev\wXwi,wj)eB 



V + 



^7 



(1) 



(2) 



2d{v,lij) 

The proof is by contradiction. 

If p ^ Cl-L{y) is closest to v ^ int CT-L(y) then, since p ^ CHiV)^ the ellipse 9p{v) must intersect 
dCT-L{y). Because p is not closer to any site in dCT-L{V) than to then 9p{v) cannot contain any 
G VF, and therefore 6p{v) must intersect some segment {wi,Wj) G >B as in Fig. 7(a) (to see 



w 



this note that the intersection between an ellipse and a line is an open line segment, and that 
the only way in which this segment can overlap WjWj but not contain either Wi,Wj is for it to be 
{Op{v) n Wj/Wj) C WjWj). 

Consider the two open half spaces and on either side of the supporting line of a 
boundary segment {wi, wj) G 0, where H^- is chosen to be on the "empty" side {H^ nCT-l{V) = 0, 



by Lem. 6.4). Clearly, {H^-,H--Jij} partition R^. Since 9p{v) must intersect some boundary 



segment, assume w.l.o.g. that it intersects {wi, wj). We show that p cannot be in either or the three 
sets in the partition, a contradiction. 

Clearly, if p G kj, then the fact that 6p{v) H kj C WjUJ] implies p G WjUJ], which is precluded by 

the fact that p is outside CH{V). 

Assume p G H^-, as in Fig. 7(b) Since not all sites are colinear, then there is ^ W that 

and therefore G H^-. Eq. [ijhas been constructed so that, for any 

IIpII > p, if m G WjWj is the closest point to p 



{■■, as in Fig. 7(b) 

is not colinear with Wi^Wj^ and ther^....^ ^ ^ 

choice of Wi^ wj^ Wk G VF, Lem. |6. 3] implies that, since 

in lij, then it is D{wk^p) < D{m^p). The convexity of D{-,p), and the fact that Op{v) Dl 
implies D(wk,p) < D(m,p) < D{v,p)^ contradicting the fact that the site closest to p is 



Lastly, if we assume that p G (Fig. 7(c)), we obtain a similar contradiction using Eq. 2 



The ellipse 9p{v) is transformed by Mp into the circle 9^ ,{v^). Since it is 9p{v) H kj C WiWj, then it 



18 



is 0'^,{v') dT-j C w'-Wj. a simple calculation reveals that this can only be true if the radius of 0'^f{v') 
satisfies 



\\p'-A\< 



However, from Eq. [2] and the fact that \\p\\ > p, it is 

\\p -v\\ > \\p-v\\ > \\p\\ - \\v\\ > ^ II ^ ^11/ V . / 



2d{v,lij] 

f Il2 lA ^(.J V ^2 



> 



|^^-^;■||74-rf (^;^^^)- 



a contradiction. 

Since we have shown that p ^ lij U H^- U H^-^ where {H^-^ H^, lij} partition R^, this creates the 
contradiction that we were after. Therefore, given the above choice of p, every p with ||p|| > p is 
closer to the boundary sites W than to any site y \ in the interior of CT-L(V). □ 

Appendix D 




Figure 8: Diagram for the proof of Lem. 2.4 



Lemma 2.4 G and P are simple (have no multi-edges or self-loops). 



Proof. By definition, G doesn't have self-loops since primal edges in P are connected sets of equidis- 
tant points to distinct sites and therefore they always connect distinct vertices. The only 
possible exception are self-loops in P connecting to itself. Because an edge of P is a connected 
set Ey^^ of points equidistant to two sites v,w^ if the edge is a self-loop of p^o then E^^^ is un- 
bounded on the two open half-spaces on either side of the supporting line ly^^ of w and therefore, 
by Lem. |6.4[ neither half-space can contain any site. Therefore all sites are contained in and 



are thus all colinear, a contradiction. Therefore there are no self-loops. 

We prove by contradiction that G has no multi-edges, and therefore, by duality, that P doesn't 
either. 



Consider two edges in the embedding G that connect the same vertices v^w (Fig. 8(a)). One 

) and consists of a point p G R{v) H R{w), a path C R{v) 



from w to p. The second edge, shown in red, has analogous 



such edge is shown in blue in Fig. 
from V to p^ and a path ^^j^ c R{i 
structure. 

From the construction of G, it is clear that the two edges must not intersect (other than meeting 
at the endpoints), and therefore their union is a simple, closed curve whose complement has, by the 
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Jordan curve theorem, an interior /. / must contain some regions corresponding to sites other than 
v^w, say zi, . . . ,Zm5 OT else the two edges would be connected and thus not counted as separate 
edges. 

Since the diagram is orphan-free, the regions inside / contain their generating sites and thus 
there are sites zi, . . . in the interior of /. We show that this is not possible, resulting in a 
contradiction. 

Consider any site z in int /, as in Fig. |8(b)[ We first split int /, as in the figure, into two 
regions A and and show that z cannot belong to either. Although A and B may not always be 
disjoint, it is easy to show that they cover int / (i.e. int I C AU B). 

A (shaded in Fig. 8(b)| ) is the set of points in / that are inside some segment vf with r G "Jvp^lvq^ 
or inside some segment wf with r G ^wp U "fwq- If a site z belongs to A then, without loss of 
generality, there is tq G ^vp such that z G vtq (the argument would be the same for 'Jvqi'ywpi'ywq) 
and so z = (1 — A)i; + Aro, A G (0, 1). Since tq G ^vp C R{v), it is D(v^ro) = minD(u^ro). However, 

z = (1 — X)v + Aro, A G (0, 1), implies that D(z^ tq) = XD(v^ tq) < D{v^ rg), a contradiction. 

B is the interior of the triangles vwp and vwq. We will only show that z ^ int vwp^ but the exact 
same argument proves that z ^ int vwq. Because p G R{v) H R{w), p is closest, and equidistant 
to v^w, and therefore the open ellipse 0p{v) cannot contain any site (or else p would belong to its 
corresponding Voronoi region instead). If some site z is z G int vwp then, since int vwp C Op{v)^ 
it is z G int vwp C Op(v), a contradiction. 

Since z ^ AU B, and int I C AU B, then it is z ^ int /, which is the contradiction that we 
were after in the first place, and therefore no two edges in G connect the same vertices. □ 



Lemma 



4.9 G has no degenerate (null area) elements. 



Proof. Because G is dual to P, to every degenerate face of G with vertices v^ w corresponds a 
primal vertex in P: a point c G R{u) H R{v) H R(w), and therefore all of v^ w are in the boundary 
of the ellipse 6c{u) (with 6c{u) = Oc{v) = 9c{w)). 

Since a line can intersect an ellipse at most at two points, three or more points cannot be 
both colinear and in d9c{u). If u,v,w are in a degenerate face, they are colinear, and in d0c{u), a 
contradiction. Therefore no faces of G are degenerate. □ 



Appendix E 



Lemma |4.1 To every boundary edge {vi^ vj] 
[B C B] 



ofG corresponds a segment in the boundary ofCT-L(V). 



Proof. By the definition of G, to every boundary edge (i;^, vj) ^ B corresponds a primal edge in P 
that is incident to the point at infinity p^o G V^. In turn, to this edge corresponds an edge in P: 
an unbounded set Ey.^y. of points closest to Vi^vj. 

Consider the two open half-planes and on either side of the supporting line kj of Vi, vj. 

it is Ey.^y.nlij C {(vi + Vj)/2}, it must be either that Ey.^y.nH:^ or Ey.^y.r\H^- 



2.1 



6.4 



they cannot both be, or else H^- DV — 



Since, by Lem. 
are unbounded. By Lem. 

therefore all sites would be in Uj (all cohnear). Assume w.l.o.g. that E^.^^^ H H^- is unbounded. 
By Lem. 



and H V = 

''J 



and 



6.4 



n V 



Eyy. n unbounded implies 
VjVj C dCl-Liy) (vi, Vj are vertices in the boundary of C1-L{V)). 



and therefore Vo 



G W and 



It only remains to show that Vi^vj are consecutive in the sequence {wi 



l,...,m). We 



prove this by contradiction. If they are not consecutive, since ViVj C dCT-LiV)^ there must be a 
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site w G VjVji 7^ Pick some p G Ey.^y., by definition closest to Vi^vj. By the convexity of 

D{-,p), it is D(w^p) < D(vi^p) = D(vj^p)^ a contradiction. 

Since 1;^, 1;^ are consecutive vertices in (wi : i = 1, . . . , m), then (vi, vj) ^ B. □ 



Lemma 



4.3 



If^{p) — {^i + with p G M^; Wi,Wj G VF; t/ien D{wi,p) = D{wj,p). 



Proof. Let m^j = (it;^ + '^j)/2. If 3:^(A) = (1 — A)it;^ + Ait;j, then the fact that 7r(p) 

that rriij = 3:^(1/2) is closest to p in the support line of Wi^wj. Therefore, ^^^^j^^'^^ (1/2) = 0, or 
equivalently: 

dDixiX) v)'^ 

2 (1/2) = (Mpit;^ + Mpit;^- - Mpp)\MpWj - MpWi) 

dX 

{MpWj - Mpp)\MpWj - Mpp) - {MpWi - Mpp)\MpWi - Mpp) 
= D{wj,pf -D{w^,pf ^0 

and thus D(wi^p) = D(wj^p). □ 



4.4 



There is p such that, for any segment {wi, wj) G every p G with \\p\\ > p whose 



Lemma 

closest point in lij is nip G WiWj is closer to V \ {wi^ wj} than to l^j. 

Proof. For each choice of {wi^ Wj) G 0, define Sij to be the set of points p G whose closest point 
in lij is rrip G WjWj. Pick some v G H C \ {wi, ^j}^ which always exists since not all sites 
are colinear. We can use Lem. 6.3, where Wi^ wj^ lij^ Sij take the role of /, 5, to conclude that 
there is a sufficiently large pij such that all p G Sij with \\p\\ > p are closer to v than to rup. By 
defining p as the maximum of pij over all (wi^wj) G 0, the lemma follows. □ 



Lemma |4.5| Every continuous function F : §^ ^ that is not onto has a fixed point. 

Proof. Assume F misses p G S^, and let 7 : §^ \ {p} be a diffeomorphism between the 

punctured sphere and the open unit disk. Since 7 o F is continuous and is compact, then the 
set C = (7 o F)(S^) C D"" is compact. 

The function g : C ^ C with g — ^ o F o 7"^ is continuous and therefore, by Brouwer's 
fixed point theorem fT3], has a fixed point x G C The fact that (7 o F o ^~^){x) — x implies 
F{'^~^{x)) = ")~^{x) and thus ^~^{x) G is a fixed point of F. □ 



Lemma 4.6 To every segment {wi,Wj) in the boundary of CHiV) corresponds a boundary edge of 



G. [B D B] 

Proof. Pick a sufficiently large p > max-^^y \\v\\ such that every p with ||p|| > p is outside CT-L(V) 
and such that Lemmas 6.5 and 4.4 hold. For any a > p, if A : C{a) C(a) is the antipodal map 
A{p) = — then, by continuity of tt and by the continuity of z/cr, the function AoV(jOti : C{a) C(a) 
is continuous. By Lemmas 4.3 and 6.5, if for some pij G C{a) it is 7v{pij) = (wi + ujj)/'^ with 
{wi^ Wj) G >B, then pij is (strictly) closest to Wi^ wj^ and therefore belongs to the primal edge E^.^^., 
which implies that {wi,Wj) G B. 

Showing that B C B now reduces to showing that for all (wi^wj) G for all a > p, there is 
Pij G C{a) such that 7v{pij) = {wi + Wj)/2. 

Assume otherwise: that there is {wi^ wj) G B such that no p G C(a) satisfies 7r(p) = (wi + Wj)/2. 
Then the function A o o tt : C{a) C{a) is not onto and therefore, by Lem. |4.5| (and using the 
fact that C{a) is isomorphic to 5^), it must have a fixed point q. 

Since {A o Ucr ^ ^){q) — Q then (z/cr o ^){q) — —Q- Since x = 7v{q) is the closest point to q in 
dCT-L(V)^ there is a segment {wk^wi) G B such that x G WkWi. Consider two open half spaces 
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and on either side of the supporting hne of Wk,wi. Since not ah sites are cohnear, we can 
choose these half spaces so that H^iHV = (p and H^iDV ^ (j). By the definition of z/cr, and recahing 
that the chosen origin of M? is in the interior int CT-L{V) of the convex huh, it is i^aix) G 
and q — —V(j{x) ^ Hj^^. To see this note that the outward-facing normal n{x) is defined so that 
X + n{x) G and so = a • n{x)/\\n{x)\\ G On the other hand, since the origin is in 

int CT-L{V)^ the fact that i^aix) E H^^ implies q — —U cr{x) G 



Since p was chosen sufficiently large for Lem. 4.4 to hold, and q G Ht^j, q is closer to some site 



V \ {wk, wi} than to wj^wi. Since v G CT-L{V)^ this contradicts the fact that x = 7r(g) G WkWi is 
the closest point to q in and thus B C B. □ 



Appendix F 

Lemma |5.4| Given a non-vanishing one-form the sum of indices of interior vertices of G is 
non-negative. 

Proof. Given non- vanishing the index of a face / is ind^^ (/) < 0. To see this, assume otherwise: 
a face with vertices vi^ . . . around it, and index one satisfies, by the definition of index and of 
n^vi < • • • < n^Vjn < n^vi (or n^vi > • • • > n^Vm > n^vi)^ a contradiction. 



Because, by Lem. 4.8, the boundary edges of G form a (not necessarily strictly) convex, simple 
polygonal chain then, given any non- vanishing all the boundary vertices have index zero, except 
for the "topmost" (argmax i^{v)) and "bottommost" (argmin i^{v)) vertices, which have index 

one. 

Since face indices are non-positive, and the sum of indices of boundary vertices is two then, by 



Lem. 5.3, the sum of indices of interior vertices must be non-negative. □ 



Lemma |5.5 // G has an edge foldover then there is a non-vanishing one-form such that 



ind^^ {v) < for some interior vertex v \ W. 

Proof. If edge foldover e = (v^ w) is a non-boundary edge, then at least one of its incident vertices, 
say V is an interior vertex v ^ V \ W. 

Consider the two faces /i, /2 incident to e, which, by definition of edge foldover, are on the same 
side of its supporting line, and the two edges ei, 62 in /i, /2 respectively, incident to v. Taking the 
half-line h from v towards w as reference, consider the (open) set Li C of directions ranging 
from h to Ci. The set L = Li H L2 is not empty since, by Lem. |4.9| /i, /2 are not degenerate, and 
therefore neither ei, 62 are parallel to h. L is also uncountable, since it is a range of the form 

L = {n G §^ : n±^h < A n±^ei > A n±^e2 > 0} 

where /i,e^ are the direction vectors of /i,e^, and n± is one of the two orthogonal directions to n, 
chosen to fit the definition. 

Because L is not empty, and it is uncountable, and because the set of edges E is finite, then 
there is always some direction n ^ L that is not orthogonal to any edge in E. We prove that the 
non- vanishing one- form is such that ind^^ {v) < 0. 

The (cyclic) sequence of oriented half-edges around v is, without loss of generality, S = 
[{v^vi); (v^w); {v^V2); . . .], and therefore the values of the one-form around v are [^^('^1) — C^('^), 
^""{w) - C^(i;), C^(i;2) - ^""{v), ...]. By the definition of n, it is ^""{vi) < ^^(1;), ^""{w) > 
^^(v)^ and ^^(^^2) < C^('^)5 and therefore the number of sign changes in the subsequence = 
[{v,vi)] {v,w)] ('^,'^2)] is four. Since the number of sign changes in the full sequence S cannot be 
less than that of its subsequence 5^, it is sc^^ {v) > 4 and therefore ind^^ {v) = 1 — sc^^ {v)/2 < 0. □ 
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Lemma 5.6 Given n ^8>^ and non-vanishing one-form , if G has an interior vertex v ^V\W 
with index ind^^(i;) = 1, then there is a face f of G that does not satisfy the empty cir cum- ellipse 
property. 

Proof. We must prove that there is a face / all of whose circumscribing ellipses contain some vertex 
in its interior. 

Consider the vertex v \ W with ind^n('^) = 1 — sc^n{v) = 1, and thus with sc^n('^) = 0. 
If [ui^U2^ . . . ^Ujn] is the cyclic sequence of vertices neighboring then sc^^i^) — implies either 
i^{ui) > C^('y), i = l,...,m, or ^'^{ui) < ^^{v), i = l,...,m. Assume the former w.l.o.g. The 
line / = {x G : n^x = ri^v}^ passing through v, strictly separates v from the convex hull of its 
neighbors. 

Consider the mesh G\ with the same structure as G but in which all the incident faces to v 
are eliminated. We show that, in G', the face count of v (the number of faces in which v lies) is at 
least one. Since / separates v from its neighbors, it also separates all the faces incident to v from 
V (except for v itself, which lies on /). Pick any direction d G 5^ with n^d < 0. The half-line h 
starting at v with direction d does not intersect any face in G that is incident to v. Since there 
is only a finite number of edges and vertices, it is always possible to choose h not to contain any 
vertex other than and not to be parallel to any edge. Since CHiV) is bounded and h isn't, there 
is some point x ^ h outside CT-L(V), whose face count must be zero. Moving from x toward h 
crosses dCT-LiV) only once (since CT-LiV) is convex), incrementing the face count to one. Because 
every interior edge is incident to exactly two faces, every subsequent edge cross (which is transversal 
because h is not parallel to any edge) modifies the face count by either zero, two, or minus two. 
Since the face count cannot be negative, and it is one at /i H dCT-L{V)^ then it must be at least one 
at V. Since G^ does not contain any face incident to this implies that there is some face / not 
incident to v such that v ^ f. 

We prove that the face / above cannot satisfy the ECE property. Since v is in. f but is not 
incident to it, and / is convex then, by Caratheodory's theorem, v can be written as a convex 
combination v = Xiui + X2U2 + X^u^^ Yli=i — ^1 K ^ (0, 1) of vertices ui^U2^u^ incident to /. 
Given an ellipse 6 circumscribing the vertices incident to /, because 9 is convex, and ui^U2^us lie 
in its boundary, then any convex combination of them with G (0, 1) must be in the interior of 0, 
and therefore / does not satisfy the ECE property. □ 



Lemma 5.7 G has no edge foldovers. 



Proof. Assume G has an edge foldover. By Lem. |5.5[ there is a non-vanishing one-form such 
that some interior vertex v G V \ W has ind^^('u) < 0. Since, by Lem. 5.4, the sum of indices of 



interior vertices is non-negative, then there must be at least one interior vertex u ^ V \W with 
positive index ind^^(i/) = 1. In that case, by Lem. |5.6| there is a face of G that does not satisfy 



the ECE property, raising a contradiction. Therefore G has no edge foldovers. □ 



Lemma 5.8 If its (topological) boundary is simple and closed, then the straight-line dual G of an 



orphan-free diagram, with vertices incident to at most three sites, is an embedded triangulation. 

Proof. Given a point p G int CT-L{V) in the interior of the convex hull of V, we show that its face 
count (the number of faces with straight edges that contain it) is one. Consider a line / passing 
through X that does not pass through any vertex of G, and is not parallel to any (straight) edge. 
It is always possible to find such a line since the set of vertices and edges is finite. Because the line 
is unbounded and CT-LiV) is bounded, there is a point x G / that is outside CT-LiV). At this point 
clearly the face count is zero. Moving from x toward / crosses the boundary of CT-L(V) (and 
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therefore, by Thm. 4.7, the boundary of G) only once, since it is a simple convex polygonal chain. 



incrementing the face count by one. At every edge crossing (which is transversal by the choice of 



line), the face count remains one since, by Lem. 5.7 there are no edge foldovers, and thus every 



non-boundary edge is incident to two faces that lie on either side of its supporting plane. Therefore 
the face count at p must be one. 

Since every point inside CT-L{V) is covered once by faces in (5, and the boundaries of G and 
CT-l{V) coincide, then (5 is a single-cover of C7-l{V). Because two straight edges that cross at a 
non-vertex always generate points with face count higher than one, then the edges of G can only 
meet at vertices, and therefore G is embedded. □ 

Appendix G 

So far we have assumed that Voronoi vertices are incident (equidistant) to no more than three sites. 
In general, however, Voronoi vertices may be incident to three or more sites, and the straight-edge 
dual G will be a polygonization, instead of a triangulation (simplicial complex). We show here that, 
even in this, more general case, the polygonization G is embedded, and can be easily triangulated 
into an embedded simplicial complex. The argument is quite simple and is given in summary. 



Theorem 5.9 If its (topological) boundary is simple and closed, then the straight-line dual G of 



an orphan-free diagram is an embedded polygonal mesh with convex faces. 

Proof. First assume that only Voronoi vertex c is incident to more than three sites {vi : i = 1 . . . m}, 
m > 3, ordered in, say, clockwise order around c (which is possible since they are equidistant to 
c and therefore lie at the boundary of the ellipse Oc{vi)). Note that, since the {vi : i = 1 . . . m} 
are equidistant to c, then any polygon (or triangle) connecting any of the {vi} satisfies the empty 
circum-ellipse property (with empty circum-ellipse Oc{vi)). The Voronoi regions {Ri : i = 1, . . . , m}, 
corresponding to those sites, are the only ones incident to c. We can order these regions in clockwise 
order around c: {i?^^ : k G K} (for some index sequence K with k ^ K satisfying A: G {1, . . . , m}), 
and therefore the polygon 11 dual to c will connect vertices {vij^ : k G K}, in that order. Now choose 
any vj in E and triangulate 11 in a fan arrangement centered at vj. As pointed out above, both 11 
and every triangle in the triangulation of 11 satisfies the empty circum-ellipse property. Therefore 
the new, triangulated G (which replaces 11 by its triangulation), satisfies all the conditions needed 



for Theorems 5.9 and |4.7| (as well as Lem. 4.9) to hold, and in particular, it is embedded. 

The claim now easily follows: any closed polygon with vertices {vi^^ : k G K}^ such that all of 
its fan-arranged triangulations centered at every one of its vertices is embedded, must be convex 
(given a polygon with non-convex vertex triangulate it in a fan centered at a vertex incident to 
w: this triangulation is not embedded). Lastly, since vertices {vi : i = 1 . . . m} are equidistant from 
c, the only polygon connecting them that is also convex is the one where vertices are arranged in 
clockwise (counter-clockwise) order around c. 

The argument has been made for a single Voronoi vertex incident to more than three sites, but 
the same argument applies to diagrams in which there are more such vertices (the key being that 
the triangulations of each dual polygon are independent of one another). □ 

Appendix H 



We show here that we may drop the bounded anisotropy condition from Thm. 6.1, at the expense 



of loosing the property of covering the convex hull of the sites; thus concluding that the dual of 
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an orphan-free anisotropic diagram is an embedded polygonal mesh with convex faces. The proof 
proceeds by reducing this case to that of Thm. |6.1[ Because of its simplicity, it is given here in 
summary. 

Assume given a continuous metric Q defined over R^, and a set V of sites forming an orphan- free 
anisotropic Voronoi diagram Vq. For now, assume that the set of Voronoi vertices is finite, and 
therefore bounded. Since there is p > such that all Voronoi vertices are inside the origin-centered 
ball B(0; p) of radius p, we can construct a new metric Q\p) = Q{p) • (1 — + / • A(p), where 



Hp) = <^ 



0, if IIpII < p 
\\p\\-p, ifp<|H<2p 

1, if \\p\\ > 2p 



Clearly, is continuous, since both Q and A are. Since it is = / outside of 5(0; 2p), and 
B{0;2p) is compact, then has bounded anisotropy (ratio of eigenvalues). We may then apply 
Thm. 



6.1 



to conclude that the dual Gg/ of Vg/ is an embedded polygonal mesh with convex faces. 
Since B{0; p) includes all Voronoi vertices of Vg, and it is = Q inside 5(0; p), then all Voronoi 
vertices in Vq are also in Vg/. Therefore, by duality, all faces of Gq are also in Gg/. Finally, since 
(5g/ is embedded with convex faces, and removing faces from a polygonal mesh preserves both 
properties, then Gq is also embedded with convex faces. 

The only assumption we have made is that the set of Voronoi vertices of Vq is finite, which can 
be justified as follows. 

First, we show that Voronoi vertices cannot be isolated, and are always incident to some Voronoi 
edge. Consider a Voronoi vertex Cijk closest and equidistant to sites Vi^vj^Vk (the reasoning for 
vertices equidistant to more sites is analogous). Since Q, and therefore D, is continuous, there 
is a (possibly small) closed ball B(cijk^e) centered at Cijk where all points are strictly closer to 
Vi^Vj^Vk than to any other site, and therefore in B{cijk^e), the Voronoi diagram of V is the same 
as that of {vi^vj^Vk}. We consider the Voronoi diagrams {vi,Vj}Q and {vj,Vk}Q of {^i^'^j} and 
{^j^^k} (restricted to B{cijk^e))^ respectively. By Corollary |2.5[ the only Voronoi edge of {vi,Vj}Q 
is connected, and likewise for {vj,Vk}Q. Since Cijk is in their intersection, then it must be incident 
to both, and therefore it is not isolated. 



By Corollary 2.5, the edges (connected sets of points equidistant to two given sites) of Vq are 
unique (thus finite since V is finite), and by the induction argument of Sec. [2| each edge connects 
two Voronoi vertices. If there are rig Voronoi edges then, since Voronoi vertices are always incident 
to some Voronoi edge, and each edge connects two vertices, there cannot be more than 2ne Voronoi 
vertices. In particular, the number of Voronoi vertices is finite. This concludes the proof. 

Finally, note that, from the way we have defined the dual (5g, and the fact that it is embedded, 
we can conclude that Voronoi vertices are unique (otherwise, multiple Voronoi vertices would result 
in coincident dual polygons, which would contradict the fact that Gq is embedded). This, together 
with Corollary |2.5| and the orphan- freedom assumption means that, as stated in Corollary |6.3| 
orphan-freedom is sufficient to ensure that all the elements (vertices, edges, faces) of the primal 
diagram are unique. 
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